Effective Sublattice Magnetization and Neel Temperature in Quantum 

Antiferromagnets 
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We present an analytic expression for the finite temperature effective sublattice magnetization which 
would be detected by inelastic neutron scattering experiments performed on a two-dimensional 
square-lattice quantum Heisenberg antiferromagnet with short range Neel order. Our expression, 
which has no adjustable parameters, is able to reproduce both the qualitative behaviour of the 
phase diagram M(T) x T and the experimental values of the Neel temperature TV for either doped 
YBa2Cu30e.i5 and stoichiometric La2Cu04 compounds. Finally, we remark that by incorporating 
frustration and 3-D effects as perturbations is sufficient to explain the deviation of the experimental 
data from our theoretical curves. 
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Two dimensional quantum antiferromagnetism has 
been a matter of great interest and subject to intense 
investigation, due to its possible relation to the nor- 
mal state properties of high-temperature superconduc- 
tors. There is by now clear experimental evidence that 
the pure high-T c superconducting cuprate compounds 
are well described by a quasi two-dimensional S = 1/2 
Heisenberg antiferromagnet on a quasi-square lattice, 
whose sites are occupied by Cu ++ magnetic ions. The 
dynamical structure factor of the 2D Heisenberg antifer- 
romagnet, calculated via the mapping of the Heisenberg 
Hamiltonian onto the 0(3) nonlinear sigma model [Q, 
was successfuly confirmed by inelastic neutron scattering 
experiments on La2Cu04 Several other microscopic 
techniques like light scattering j3j, muon spin relaxation 
Q and thermal neutron scattering ||, have also been 
used to probe the magnetic correlations in these materials 
and confirmed the quasi 2D Heisenberg antiferromagnet 
hypotesis. 

A common feature among almost all superconducting 
cuprate compounds is the existence of a Neel ordered 
moment in the low temperature, underdoped regime. As 
the temperature is increased, or the sample doped, an- 
tiferromagnctic order is destroyed, leading to new forms 
of spin order |q|. According to spin- wave theory, for a 
d-dimensional hypercubic lattice, Neel order is possible 
at T = for d > 2. However, despite the widespread suc- 
cess of spin- wave theory, there remain a number of issues 
that defy the description of the superconducting cuprate 
compounds within this approach. For example, the de- 
velopment of a sublattice magnetization is known to be 
suppressed in the two-dimensional Heisenberg antiferro- 
magnet for any nonzero temperature. True long range or- 
der, as a genuine three-dimensional phenomenon, would 
only be achieved by considering the interlayer coupling 
Jj_ ~ 10~ 5 Jii not only as a perturbation. 



It is the purpose of this work to show that the experi- 
mental data for the sublattice magnetization of La2Cu04 
and YBa2Cu3C>6.i5 || can in fact be described still 
in the context of a two-dimensional square-lattice quan- 
tum Heisenberg antiferromagnet at finite temperatures, 
as far as inelastic neutron scattering experiments are con- 
cerned. Our starting point is the observation that the 
nature of the spin correlations in the renormalized clas- 
sical regime is consistent with one of the three possibil- 
ities of fig. [I], according to the observation wave vec- 
tor \k\, or frequency ui In this sense, any possible 
neutron scattering experiment, with high enough energy 
transfers, performed on a true two-dimensional system, 
would actually measure a nonvanishing effective sublat- 
tice magnetization, since one would be probing the dy- 
namics of spin correlations in the intermediate Goldstone 
region. Inelastic neutron scattering experiments probe a 
microscopic, short wavelength physics to which we can 
associate an effective Neel moment. As it will become 
clear, the behaviour of this effective moment can be de- 
scribed by an effective field theory for the low frequecy, 
long wavelength fluctuations of the spin fields about a 
state with short range Neel order. We will then be able 
to speak about a finite temperature phase transition in 
the 2D system, associated to the colapse of the Goldstone 
region in fig. [I]. 
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FIG. 1. Properties of the 2D quantum Heisenberg anti- 
ferromagnet as a function of the frequency u>. £ is the actual 
correlation length while £j is the Josephson correlation length 
related to the spin-stiffness by £j = hc/p s - 
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The two-dimensional square-lattice quantum Heisen- 
berg antiferromagnet has a well known continuum limit 
given in terms of the 2 + 1 dimensional 0(3) quantum 
nonlinear sigma model JTof . The later, on the other hand, 
is defined by the partition function 



where the action 



T>i%i S(nf 



r 8% 



1) exp(-I(n,)), (1) 



(Vn ; ) 2 + 4(£U0 2 1 (2) 



describes the long-wavelength fluctuations of the stag- 
gered components of the spin-field rii — (er, if ) , I = 
1, . . . , N = 3. The fixed length constraint is understood. 
In the above expression, po is the spin stiffness, cq is the 
spin-wave velocity, (3 — (fe^T) -1 and all quantities with 
a subscript represent bare quantities. 

We shall work in the natural units ks = h = c = 1, 
with c being the renormalized spin wave velocity. |^] Also, 
further analysis will be simply expressed in terms of the 
coupling constant go = N/po, which has the units of in- 
verse length. With this notation and choosing the stag- 
gered magnetization to be along the a field direction, we 
can integrate over the remaining N — 1 spin- wave degrees 
of freedom if and study the behavior of the partition func- 
tion (0) in the large N limit. As usual, N is taken to be 
large enough while go is kept fixed. This means that we 
have to choose po ~ N. 

For large N, the partition function (jl|) is dominated by 
the stationary configurations of the magnetization, (er) , 
and of the Lagrange multiplier field, i(A) = m 2 , intro- 
duced in order to ensure the averaged fixed length con- 
straint. These, on the other hand, can be determined 
from the stationarity conditions 
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where a cutoff A was introduced to make the momentum 
integral ultraviolet finite. 

From the above set of equations we see that the anti- 
ferromagnetic system could in principle be found in two 
distinct phases. If (er) ^ then m = and the system 
would be in the Goldstone phase with a nonvanishing net 
sublattice magnetization. In this case the ground state 
would exhibit true long range Neel order. If m ^ on 
the other hand, then (er) = and Neel order is absent. 
There are no gapless excitations in the spectrum of the 



*For large N the spin wave velocity does not renormalize 
and c = c. 



finite temperature system. It is a well known fact that 
for the 2 + 1 dimensional O(N) invariant nonlinear a 
model at T > 0, the only possible physical situation is 
the second one, due to severe infrared divergencies in the 
second saddle-point equation (0). As a consequence, the 
value of m is pushed from zero to a finite value making 
the sublattice magnetization (er) to vanish, in agreement 
with the Coleman-Mermin- Wagner theorem. 

We can compute the value of the O(N) invariant mass 
to in a closed form by subtracting the linear divergence 
in the second saddle-point equation (|3|) as 
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where g c = 47r/A is the bulk critical coupling and 
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with S =1/2 and a being the lattice spacing pffl . Now, 
after momentum integration and frequency sum we arrive 
at 



2 . f e -PPs/(2N) 

£ = m((3) = — arcsinh 

p V 2 
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which is nonvanishing for T > 0, thus indicating that a 
Neel phase can only occur at T = 0. The conclusion is 
that even at the smallest temperature there is a gap in the 
spin- wave spectrum and a finite correlation length which 
measures the size of clusters in which there is short range 
Neel order. The above expression for £ has been obtained 
by Chakravarty et al. |Q and its zero temperature limit 
successfully confirmed by quasi-elastic neutron scatter- 
ing experiments on La2CuC>4 jnj, for, however, temper- 
atures approaching TV from above, T — > Tjy , where this 
compound is known to exhibit a true two-dimensional 
behaviour. 

Let us now consider inelastic neutron scattering ex- 
periments, performed on a true 2D system, with energy 
transfers AE = fiuj such that the corresponding wave- 
length, A = 1/fruj, satisfies £j < A < £. Typical time 
scales in such experiments are t\ = A/c, consequently 
much smaller than the relaxation time r = £/c at which 
the 2D system disorders. For such experiments, spins 
would look like as if they were frozen and a nonvanish- 
ing effective sublattice magnetization would be measured. 
Since at low temperatures £ is much larger than £j, the 
three regions of fig. [I] are well separated. In the large in- 
termediate region, probed by our experiment, the system 
behaves as if it had true long range antiferromagnetic or- 
der and dynamic scaling hypotesis is justified jj. We are 
then allowed to apply a hydrodynamic picture for the low 
frequency, long- wavelength fluctuations of the spin-fields 
77; , in which its short- wavelength fluctuations follow adi- 
abatically the fluctuations of the disordered background 
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whose typical wavelenth is the scale of disorder, the cor- 
relation length. The effective field theory to describe the 
spin correlations in this intermediate Goldstone region is 
obtained by functionally integrating the Fourier compo- 
nents of the fields in (Q) with frequency inside momen- 
tum shells k < \k\ < A. The resulting partition function 
is such that, for large N, the leading contribution now 
comes from the scale dependent stationary configurations 
(a) K and i{A) re = m 2 ., solutions of the new set of saddle- 
point equations 
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Differently from the previous case, now we can in fact 
find the system in two different phases (regimes): ordered 
(asymptotically free) or disordered (strongly coupled); 
depending on the size of £ K = 1/k relative to £: smaller 
(high energies) or larger (low energies). In the ordered 
phase, £ K <C £, m K — is the solution that minimizes the 
free energy and the 2D system is then characterized by a 
nonvanishing effective sublattice magnetization (cr) K =/= 0, 
a divergent effective correlation length £ e f / = 1 /m K = oo 
and gapless excitations in the spectrum. 

The effective sublattice magnetization can be exactly 
computed from the second saddle point equation . Us- 
ing the renormalization scheeme defined by (Q), we ob- 
tain, after momentum integration and frequency sum, the 
expression 



2nN 



^ + J-ln(2sinh(/3 K /2)), (8) 
4ttN 2np 



which depends on the energy scale k and on the temper- 
ature. 

Some comments are in order. The running spin stiff- 
ness (||) decreases (g(K,(3) = N/p 3 (K,f3) increases) as 
£ K — > £, for a given temperature. This is a consequence 
of the fact that we are coarsing over degrees of freedom 
which actually feels the finite size of the clusters with 
short range Neel order. Furthermore, for £ K > £ we would 
be coarsing over degrees of freedom outside these clus- 
ters, leading to the disordered (strongly coupled) phase. 
Lowering the scale k is also equivalent to waiting longer 
for a response, and for t k > r we would be waiting long 
enough for the system to disorder. Here, conversely, in 
order to obtain a finite temperature phase transition in 
the 2D system, we will rather fix the scale k and study 
the behaviour of the effective spin stiffness (||) with the 
running parameter being the temperature. We must fine 
tune, and hold fixed, the energy transfers in our exper- 
iment so that our 2D system is able to reproduce the 
observed 3D behaviour in real materials. For this it suf- 
fices to impose the boundary condition 



with p s being the bulk spin stiffness of the real system. 
From (^|) we conclude that k = p g /N, which is exactly 
the inverse Josephson correlation length, k = ^J 1 . This 
should not be surprising since the spin stiffness is itself 
a microscopic, short wavelength quantity defined at the 
Josephson scale. Notice also that k = 23 meV for the case 
of La2CuC>4 , which is actually consistent with the energy 
transfers commonly used in this kind of experiment [0 . 
Now, inserting (|^) in (||), the expression for the finite 
temperature effective sublattice magnetization, M(T) = 
p s (n = p s /N,T), becomes 



M(T) = f 
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with Mq — p s . The sublattice magnetization M(T) van- 
ishes at a Neel temperature Tv given by 

Mn 



T ^MnT ^ 

For temperatures above Tv, we would be coarsing over 
degrees of freedom outside the shrinked clusters of size £, 
leading again to the disordered (strongly coupled) phase. 

Let us now show that the above analysis can in fact 
be used to describe the experimental data for different 
cuprate compounds. Take for example the data obtained 
for La 2 Cu0 4 and YBa 2 Cu 3 6 .i5 @- For these com- 
pounds, we find agreement between our predictions and 
the observed Neel temperatures, within almost 10%, al- 
ready at the leading order, as can be seen from table 
[| More important, we have obtained a good qualitative 
agreement between the phase diagram M(T) x T and ex- 
periment, for the whole range of temperatures from to 
Tjv (see dotted lines in figs. || and ||). 

In order to have a flavor on how our results can be 
improved, let us mention that already at the next-to- 
leading order in the 1/N expansion we will have a non- 
trivial renormalization of the spin-wave velocity due to 
the self interaction between spin-waves |]. This should 
lower the value of c and, if we take for example a lowering 
of about 10%, we obtain the behaviour described by the 
solid curves in figs. |^ and [|. The spin wave velocity will 
also be renormalized by dynamic scaling, but we assume 
that at the shortest distances, that is fiui 3> this 
effect can be neglected when compared to the effects of 
the self interactions. For this reason, it should not lead 
to a further damping of the spin waves. 





c (eV A/h) 


p s (meV) 


Tjv (K) 


T e N xp (K) 


YBazCuaOe.is 


1.00 ±0.05 
0.90 ± 0.05 


81 ±4 
73 ±4 


450 ± 20 
410 ± 20 


410 ±3 


La 2 Cu04 


0.85 ±0.03 
0.75 ±0.03 


68 ±2 
60 ±2 


380 ± 10 
330 ± 10 


325 ±5 



p s (K, T = 0) = p s 



(9) 



TABLE I. To compute p s we have made used of formula 
(^) with a ~ 3.8 A and S = 1/2, while to compute Tn we 
have used formula ( |TT| ) with iV = 3. The experimental values 
of c were taken from ref. fl2| . 



3 



Notice now that, with respect to the solid curves, the 
experimental points can be separated into two different 
sets. For T < Tjy/2, we find all points below the solid 
curves while, for T > TV/2, we find, instead, the points 
all above our theoretical prediction. This is consistent 
with a picture in which strong frustration induced quan- 
tum fluctuations, due for example to a nonzero next- 
nearest-neighbour coupling, are dominant at low tem- 
peratures and suppressed at higher T, where the effects 
due to a sizable J± begin to be felt. Notice also that in 
the case of YBa2Cu3C>6.i5 the points deviate even more 
from the solid curve, for T > TV/2, than in the case of 
La2CuC>4 . We attribute this to the bilayer structure of 
YBa2Cu30g.i5 , which causes a further increase in the 
sublattice magnetization. As we approach T/v from be- 
low, both systems behave effectively as true 2D Heisen- 
berg antiferromagnets with nearest-neighbour coupling, 
as shown by the experimental data. From the above 
discussion we conclude that by incorporating frustration 
and 3-D effects as perturbations, with properly temper- 
ature renormalized coeficients, might be sufficient to ac- 
count for the deviation of the data from our theoretical 
predictions. We are presently investigating this possibil- 
ity. 
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FIG. 2. For the dotted line we have used he = 0.85 eV A 
while for the solid line we assumed he = 0.75 eV A. Experi- 
mental data from ref. Jjj]. 

As a final remark, let us show that our treatment is 
consistent with experiment also in the disordered phase. 
Above T/v, where the effective sublattice magnetization 
vanishes, we must consider the second possible solution 
for the set of saddle-point equations (0), namely (a) K = 
and m K ^ 0. If we then solve the self-consistent equation 



gap we end up with 
4 



4 2 / e -/W(2iV) 

mi = -^-arcsinh 

K R2 



for k = p s /N. It is straightforward to see that the 
effective correlation length £ e // = l/m K diverges as 
T — * p s /(iVln2), or in other words, as we approach T/v 
from above. This is consistent with the data for La2Cu04 
from pd| ], and for the correct temperature limit. 
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FIG. 3. For the dotted line we have used he = 1.00 eV A 
while for the solid line we assumed he = 0.90 eV A. Experi- 
mental data from ref. |s|. 
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